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The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation 
mechanisms of high-energy emission from these systems. The high statistics allow for careful modeling of the light curve features 
as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA I pulsars with simulated 
high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models, within the 
vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore 
the phase space of the magnetic inclination angle, viewing angle, maximum emission radius, and gap width. We also used the 
measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius, under the assumptions 
that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of 
curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from 
the vacuum field more closely match the observed light curves and muitiwavelength constraints, and that the calculated parallel 
electric field can place additional constraints on the emission geometry. 
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1. INTRODUCTION 

The pulsar emission mechanism is not well understood. 
Magnetospheric particle acceleration is likely responsible 
for the observed emission, but the emission geometry is 
unknown. One can gain some insight by comparing light 
curves derived from geometrical emission models with ob- 
served pulsar light curves. Observations of pulsars by 
the Fermi Gamma-ray Space Telescope Large Area Tele- 
scope (LAT) [4] have shown that the high-energy emission 
likely originates in the outer magnetosphere [1]. We sim- 
ulated high-energy light curves from geometrical versions 
of two standard high-altitude emission models, the outer 
gap (OG) [13] and slot gap (SG) [11], and from a third, 
modified SG model with azimuthal asymmetry in emis- 
sivity due to a naturally occurring offset dipole (aSG) [8], 
These models were considered within two field geome- 
tries, the vacuum retarded dipole (VRD) and force-free 
(FF) [6] fields. We compared the resulting light curves 
with the LAT light curves of the Vela pulsar and PSR 
J0007+7303, the CTA 1 pulsar, to constrain the systems’ 
geometries, and calculated the model-dependent magni- 
tude of the accelerating electric field fit in the VRD field. 


2. LIGHT CURVE MODELING 

To model the LAT light curves, we first simulated pul- 
sar light curves from geometrical representations of the 
SG, aSG, and OG emission zones within the VRD and FF 
fields, following the simulation method of [7], The i? field 
defined in the observer’s frame is transformed to the co- 
rotating frame (CF) [5 ] and photons are emitted tangent to 
B in the CF prior to calculation of the aberration. We as- 
sume constant emissivity along the field lines in the CF. 
The azimuthal asymmetry in the polar cap (PC) angle for 


the aSG model is calculated for each inclination angle a 
as in [8] (see also [9]). For a given a , gap width w (in 
units of open volume coordinates r ovc , as in [7]), and max- 
imum emission altitude r, the code outputs the dimension- 
less emission intensity and the minimum and maximum 
radii of curvature, p min and p max . emission radii r mm and 
r max , and local field magnitude |5| m in and |S| max , at all ob- 
server angles < and rotation phases <p. 

We simulated light curves for a fiducial rotation period 
of 0.1 s on a 4-dimensional grid of a, C, w, and r. Our 
simulation resolutions are 1 ° in a for the VRD field and 
1 5° for the FF magnetosphere; 1 ° in 0.01 r ovc in 0 r ovc < 
w < 0.3 r ovc ; and 0.1 Ri c in 0.7 R\ Q < r < 2.0 R k , where 
R k = c/Q is the light cylinder radius. Emission is allowed 
out to a cylindrical radius r cy i = 0.98 R k for the OG model 
and r cy i = 0.95 R\ c for the SG models. 

LAT light curves were constructed from photons in 
an angular radius 0 < max[1.6 - 31og m (£), 1 .3] from 
the pulsar [2], The PSR J0007+7303 light curve has 
32 fixed-width bins, and was taken from [3], The Vela 
light curve has 140 fixed-count bins of ~ 3000 pho- 
tons each. The background of PSR J0007+7303, 195 
counts/bin, was found using the Fermi tool gtsrcprob as 
in [3], The emission in the off-peak above the background 
level was assumed to be magnetospheric in origin (but 
see [3] for details). The Vela background was found to 
be 204 counts/bin using the off-peak (phases 0.8-1 , where 
no magnetospheric emission was detected in our spectral 
fits) counts in the energy-dependent PSF of the pulsar. 

We used a Markov Chain Monte Carlo (MCMC) maxi- 
mum likelihood routine [14] to search the parameter space 
for the combination of (a, <T, w, r, A <f>) that best repro- 
duced the LAT light curves. The fifth parameter, A <p, is the 
amount by which a model light curve must shift in phase 
in order to best match the LAT light curve. The MCMC 
begins at a random point in parameter space, calculates the 
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Figure 1 : Light curve modeling results for Vela and PSR J0007+7303. The absolute best fit parameters are given next to each model; 
uncertainties in individual parameters vary, but are typically of order - 10%. In all panels, the light curve is shown in black. Blue 
curves show the best fit OG light curve, red the best fit SG light curve, and green the best aSG light curve. Vertical dot-dashed lines 
show phase zero of the model light curve of corresponding color: the phase of the line is the same as the listed A <f>, and is the phase of 
emission coming from the closest magnetic pole. The horizontal purple dashed lines show the background count level. Panel (a) 
shows results for the Vela pulsar with the VRD field, (b) Vela with the FF field, (c) PSR J0007+7303 with the VRD field, and (d) PSR 
J0007+7303 with the FF field. 


likelihood L , and moves in a random direction to a new 
point in space to calculate L again. If L/Lf ormer > 1 or 
> rand[0, 1), L is saved in a chain; the routine runs until the 
chain contains the user-specified number of steps. Many 
chains were run to explore the whole parameter space. We 
used Wilks’ theorem, A In L = -Ay 2 / 2, to calculate L at 
each point in parameter space. To perform a fit, we sub- 
tracted the background level from the LAT light curve, re- 
binned die model to match the data, and normalized the 
model to the total counts in the LAT light curve. For each 
fit in the vacuum case, we ran 200 chains with 20 steps 
each to adequately sample the parameter space. For the FF 
fits, we ran 20 chains of 20 steps each for each a. Multi- 
wavelength constraints on a and £ were considered after 
fitting. Our fit results are shown in Figure 1. Each plot 
shows three instances of the LAT light curve with the best 
OG, SG, and aSG light curve superposed. The zero phases 
of each model, corresponding to the nearest magnetic pole, 
are given by the vertical lines, and the background levels 
with horizontal lines. 

We find that for the (a) VRD and (b) FF fields, the OG 
model (blue) statistically fits the Vela light curve better 
than either SG model. The SG (red) produces too much 
off-peak emission, leading to high y 2 values. The aSG 


(green) reduces the background significantly, leading to a 
much better fit. The aSG also qualitatively fits the peak 
emission well, as its main peaks are the correct approx- 
imate height and an inner peak is present. All three fits 
have C close to the value determined from the X-ray torus 
geometry [10], ( ~ 64°. Ad> is consistent with the observed 
phase lag between the radio and y-ray peaks for the VRD 
models; however, for the FF case, A</> is too large. 

For the pulsar in CTA 1 , the SG models fit much better 
than the OG in both field geometries. There is no con- 
straint on A <p due to the lack of a radio detection. The 
VRD geometry in (c.) produces much better fits than the 
FF in (d); there is little difference between the SG and aSG 
due to the small a. A larger FF a leads to a larger PC offset, 
which lowers the first peak. The large difference between 
a and C is consistent with the pulsar being radio-quiet due 
to geometry-the radio beam would not cross our line of 
sight for such a large |or - £|. 


3. CALCULATION OF £j. 

The y-ray spectra of pulsars are well fit by an exponen- 
tially cut-off power law. 
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where N 0 is the differential flux, Eo the energy scale, T the 
power law index, and E c the cutoff energy; b = 1 results 
in a simple exponential cutoff, while b < 1 gives a sub- 
exponential cutoff and b > la super-exponential cutoff. 
For phase averaged spectra, b < 1 due to blending of E c 
as it varies with phase (e.g. [2]), while b is consistent with 
(and is fixed to) 1 in individual phase bins. 

In current models of pulsar emission, at energies above 
~ lOOMeV the emission is dominated by curvature radia- 
tion. Particles reach the radiation reaction limit at Lorentz 
factors yc,R ~ 10'. In this limit, the curvature radiation 
cutoff energy Ecu is related to E\\ and p c by 


Eqs. - 2 p Tcr — 0-32/l c 



(2) 
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Assuming all emission with E > lOOMeV is due to 
pure curvature radiation and that the VRD is the true B 
field structure, we calculated E\\ in each light curve phase 
bin. We used the simulated minimum radii of curvature 
(Pc = Pmin) from the best fit VRD light curves of §3 and 
the measured cutoff energies (£cr = E c ) in each phase bin. 
The cutoff energies are given in [3] for PSR J0007+7303. 
We updated the Vela 0. 1-1 00 GeV phase resolved spectral 
results with 30 months of LAT data, following the method 
of [2] with 3000 pulsed counts per bin, and used our mea- 
sured cutoff energies for the Vela E\\ calculation. As an 
example. Figure 2 shows the measured E c , the simulated 
Pmin* r min , and |5| max , and the calculated £j| for the Vela 
pulsar peak emission (phases 0 < <t> < 0.8), using the best 
fit parameters from the VRD SG geometry. 

We explore how the parallel electric field varies with 
emission altitude. We have calculated E\\ in each phase 
bin for the best fit vacuum OG and SG model parameters. 
Because the value of |B|niax corresponds to r m j n , we have 
plotted E\\ and £||/|,B|max with minimum emission radius 
for the OG and SG models in Figure 3. 

For both emission models, Vela has overall a ~ constant 
or gradually varying £j| with altitude (panel (a) of Fig- 
ure 3), which is expected (e.g. [2]). In panel ( b ), the value 
of £ji is compared with |S| max . As expected for an E field 
induced by B, the ratio £i|/|,B| m ax < I for all r m i„ out to 
the light cylinder radius (near and beyond 1 R ic , the vector 
components of B are less certain and are not included in 
this calculation). 

For PSR J0007+7303, the magnitude of £j, with r min is 
consistent with a constant (panel (c)), and its values are 
similar to those calculated for Vela. Note that the geomet- 
rical parameters obtained from the light curve fits are very' 
different from those of Vela, leading to a different range of 
r min for PSR J0007+7303. The ratio £i|/|fi max < 1 for the 
best fit SG model, but it is > 1 for the best fit parameters of 


Figure 2: Example of simulated and measured values described 
in §3. (a) The LAT light curve of Vela (purple) and measured E c 
with phase (yellow diamonds), (b) Simulated p mm (solid blue) 
and r min (dashed green) for the best VRD SG fit parameters, (c) 
Simulated |B| ma * (solid cyan) and calculated £j| (red squares) for 
the best VRD SG fit parameters. 

the OG model. There are instances in, for example, a non- 
ideal magnetosphere [10] where E\\ may be > 'B\. In the 
case of the vacuum field, in which there are no currents, 
the only source of E\\ is induction by B, and therefore £j| 
cannot be larger than |5| max . In this particular case, then, 
we find that the combination of the OG model and VRD 
field we have used does not approximate the physical envi- 
ronment of the pulsar magnetosphere and/or tire geometry 
of the emission zone. 


4. CONCLUSIONS 

We have evaluated the geometries of the slot gap and 
outer gap emission models, and the vacuum retarded 
dipole and force-free magnetic field solutions, by compar- 
ing the simulated light curves with the LAT light curves 
of Vela and PSR J0007+7303 and finding the geometri- 
cal model parameters that best reproduce the data in each 
case. In general, the OG has no off-peak emission (this 
is largely responsible for the OG fitting Vela the best and 
PSR J0007+7303 the worst), while the SG models do a 
better job of reproducing wing emission but over-predict 
the off-peak emission. Introducing azimuthal asymmetry 
of the PC angle in the aSG model leads to a reduction in 
the off-peak emission level, improving upon the SG model 
light curve fits within the VRD field. In the FF field, how- 
ever, aSG light curves tend to have a much reduced first 
peak, leading to significantly worse light curve fits. 
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Figure 3: £[ (a) and the ratio of £i|/|B| max (b) for the OG (blue squares) and SG (red diamonds) models that best fit the Vela LAT light 
curve, (c) Same as (a) for PSR J0007+7303. (d) Same as (b) for PSR J0007+7303. 


The first y-ray peak occurs at later phase in FF light 
curves, so values of A <b are larger in the FF magnetosphere 
than in the VRD field. Physically, A <,6 cannot be larger 
than the phase lag between the radio and first y-ray peak 
unless the radio beam model is highly contrived. The FF 
field requires a A 4> larger than this phase lag for Vela. This 
suggests that the true pulsar magnetosphere may be signif- 
icantly different from the FF magnetosphere. 

For Vela, all fits within the VRD field have f close to 
the expected value of 64°, and all have reasonable values 
of \a - g\ such that the radio emission is observable. Inter- 
estingly, the aSG model has C closest to 64°, and while its 
X 2 is poor, it is an improvement over the SG and qualita- 
tively reproduces well the major features (two main peaks 
and inner peak) of the pulsed emission. The FF fits also 
get close to the correct £ value, but only the OG has an 
acceptable la - ( |, while the best fit SG models are consis- 
tent with a radio-quiet pulsar. Both field structures lead to 
large |a - g\ for PSR J0007+7303, consistent with the lack 
of detected radio emission. 

We calculated the model-dependent E\\ for the OG and 
SG geometries, and found that it is constant or slowly 
varying with emission radius. Comparing £y with | Blmax 
leads to the interesting result that for PSR J0007+7303, 
£jl > |.B! ma x for die OG parameters that best fit the LAT 
light curve. This is not consistent with the VRD where 
|£j < | B\, and dius disfavors the OG model in the VRD 
geometry at the location in parameter space where the best 
light curve fit is found. The calculation of E\\ can therefore 
be used as a diagnostic of the model magnetosphere and 
emission geometry, in addition to the light curve fitting. 

We cannot rule in favor of the SG or OG from these light 
curve fits. As the models are purely geometrical, we expect 
to reproduce only dominant light curve features, and our 
statistical fits are poor. The vacuum field produces better 
light curve fits than the force-free field. However, we note 
that the resolution in a of our FF models is much lower 


than in the VRD, so further modeling with higher resolu- 
tion is needed to confirm this result. We have demonstrated 
that light curve modeling leads to constraints on the geom- 
etry of individual systems, and that the phase lag is an im- 
portant diagnostic in comparing magnetic field structures. 
The model-dependent calculation of £j|, and comparison 
with the model |£j max , can additionally be used to constrain 
both the magnetosphere and emission geometry. 
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